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The dynamics of photonic wavepacket in the effective oscillator potential is studied. The oscillator 
potential is constructed on a base of one dimensional photonic crystal with a period of unit cell 
adiabatically varied in space. The structure has a locally equidistant discrete spectrum. This leads 
to an echo effect, i.e. the periodical reconstruction of the packet shape. The effect can be observed 
in a nonlinear response of the system. Numerical estimations for porous-silicon based structures are 
presented for femtosecond ThSapphire laser pump. 



One of the strongest driving forces of the modern op- 
tics is a projection of the solid-state physics concepts on a 
photonic language. Except the fundamental interest, this 
is related to so-called photonics: a construction of optical 
based elements of information processing devices. Mod- 
ern technology provides an opportunity to construct a 
photonic devices of a sunmicron and micron sizes. In par- 
ticular it is worth mentioning photonic crystals 1 ^, struc- 
tures with a periodic modulation of optical properties. 
They act on a photon the same way as a crystal lattice 
potential on electrons. Beauty optical analogs of several 
solid-state phenomena have been observed in systems, 
based on photonic crystals. One could mention Bloch 
oscillations 3,4 , optical analog of Franz-Keldysh effecl£i£, 
fabrication of optical molecule^, etc. Optical field local- 
ization in the photonic crystal based structures provides 
numerous nonlinear optical phenomena 8 ' 9 . 



A description of experimental results in this field com- 
monly uses spectral (plane- wave) representation. At the 
same time, modern femtosecond lasers provide pulses of 
a small longitude. For example 100 fs impulse has a 
30/im length in vacuum and is even shorter in a mate- 
rial media. This is comparable to the size of a photonic 
crystal structure. By taking into account a finite length 
of the photonic wavepacket, one can predict a new range 
of optical effects, similar to the ones, based on electron 
wavepacket spatial localization. In this paper we discuss 
possible observation of the "optical echo" effect, i.e. pe- 
riodical reconstruction of a shape of the optical pulse in 
a specially designed photonic crystal structure. 



The simplest echo effect can be realized for a quan- 
tum particle in a harmonic oscillator potential. Let the 
particle being localized at the initial moment near cer- 
tain spatial point xq- It can be characterized by the 
wavefunction ^(x), for example a Gaussian packet can 

be considered ^o(x) oc exp (— + ikx^j . After the 

initial moment a space spreading of the pulse takes the 
place because of the dispersion. However, after the time 
equal to the period of the oscillator T, the wavepacket 
will reconstruct it's shape back. Indeed, all the eigen- 
functions of oscillator evolve with multiple frequencies 
of the oscillator, so that the wavefunction of the packet 



^>(x,t) repeats itself with a period T: 
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One can see from (flj that the crucial circumstance here 
is that eigenlevels of the harmonic oscillator are equidis- 
tant in frequency domain. On the other hand energy 
levels of almost any pendulum are equally spaced in en- 
ergy in the quasiclassical limits. In this case energy split 
in in neighboring levels is just where T(E) is a pe- 

riod of oscillation at given energy. This means the echo 
effect can be observed for packets composed of locally 
equidistant states of the discrete spectrum of almost ar- 
bitrary potential. Echo effects can also be observed in 
more sophisticated cases, for example in the disordered 
structures 1 ^. The only requirement is the local equidis- 
tance of the spectrum. 

Our goal is to construct an effective oscillator potential 
for optical pulses. This requires a modifying of dispersion 
relation with respect to spatial coordinate. Fabricating a 
structure with a refractive index being changed smoothly 
in a wide range is a very complicated task. However, 
the structures based on one-dimensional photonic crys- 
tals can be used. In order to form a structure with a 
discrete spectrum (i.e. localized eigenstates), we propose 
using a photonic crystal with a unit cell period adiabat- 
ically varied in space. Let us imagine that for the center 
of the structure the carrier frequency of the wavepacket 
lies just above the photonic bandgap. Consider the struc- 
ture with a period decreasing to the spatial edges. In this 
case the photonic bandgap of a crystal is higher at the 
crystal edges, than in the geometric center. This means 
the wavepacket is reflected from the Bragg mirrors near 
the structure borders, therefore it becomes localized in 
the whole crystal. 

The one-dimensional photonic crystals based on porous 
silicon are quite widespread in modern photonic technol- 
ogy. The procedure for their preparation is quite sim- 
ple. Such photonic crystals constitute of repeating pair 
of layers of n ~ 1.5 and n! ~ 2.2 refraction indexes. The 
optical lengths of both layers in each repeating pair of 
the crystal are usually taken equal, their material dis- 
persion can be examined as linear. We will address to 
consider such type of structures. The adiabatic variation 
of the period of such crystal is also feasible technologi- 
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cal task. However in order to obtain sufficient number of 
equidistant localized levels (around 10) the whole num- 
ber of pairs must be taken large enough. In our model 
M = 100 pairs of layers are considered. 

In the effective local linear susceptibility approxima- 
tion all layers can be characterized by two values: d m 
- the length of the layer and n m - it's refraction index. 
The held strength of eigenmode of each layer is given by 
expression: 
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where A m and B m are unknown amplitudes, which are 
determined by frequency of eigenmode w. 

There are different ways to calculate the spectrum of 
the structure. One possibility would be a direct search 



of a solution of eigenvalue problem. The boundary con- 
ditions on the edges of the structure and on the joints 
of the layers produce a set of equations for A m , B m for 
each layer. The solvability of this set requires finding a 
determinant of a 2M x 2M matrix with some elements 
equal to zero. This would give an equation for w. The 
numerical solving of this equation produces a spectrum 
of whole system. 

Another possible way to find a spectrum of the sys- 
tem is to treat it as an effective adiabatically varied peri- 
odic potential. This method is quite similar to employing 
Vloke or Bloch's theorem. It fits the paradigm of apply- 
ing condensed matter physics ideas on the field of optical 
problems. Let's zoom on pair of two layers : (d,n) and 
(D — d,n'). Their overall spatial length is D. Each eigen- 
mode of the pair j can be presented the same way as ([2]) : 
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<d(x) is the Heavyside function. The nil of the coordi- 
nates is placed at the beginning of the first layer. Simi- 
larly to Vloke's theorem we propose that: 



£(x + D) = £{x)e i 
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Let T = k = f , and xo = kdn = k(D - d)n' . Us- 
ing the boundary conditions on the layer connections the 
value of Bloch's phase \ can be found^: 
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That is equivalent to 
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Actually, The Vloke theorem itself is not used, we just 
use the corresponding notation. In fact, the phase x is 



determined only by the parameters of local pair of layers, 
i.e. it is solely connected with the pair number. The 
expression for the value of the electric field in the pair 
number j can now be rewritten as: 



£{x) = £o(x)e 
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£q (x) is an expression for the spatial part of electric field 
in the first pair of layers. For example in odd layers (??) 
will appear as 

£,{x) = (A (k,d y kx + B (k,d )e-* kx y^ (7) 
The same equations can be derived for complex conju- 
gate £* (x) . We will consider full field strength as a sum 
of the £j{x) and £!(x) which is in fact 2$l£j(x). In this 
case J7]) becomes 
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£(x) = (A (fc,do)e Jfcx +B (k 1 d Q )e-* kx y^ + (A*(k, d a )e^ + B*(k, d o y kx )e~^ (8) 
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In order to get a spectrum of the eigenmodes of this sys- tions on the spatial edges of the system are 
tern it is virtually placed in an opaque resonator. Let L 
be a whole length of the system. The boundary condi- 



£(0) = 0,S(L) = 
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FIG. 1: The part of the spectrum of ordinary photonic crys- 
tal and a photonic crystal oscillator with 10% period width 
variation starting from the upper edge of photonic bandgap 
spectrum plotted in arbitrary K = w, c = 1 versus the num- 
ber of the level. Both structures consists of 100 pairs of layers. 
The important localized levels are highlighted. 
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FIG. 2: The difference between energy levels versus the num- 
ber of the level plotted in the same region for same structures, 
as in Fig. 1. In the frequency interval of the laser signal 
the photonic crystal with a period variation produces "more 
equidistant" energy levels than an ordinary crystal. 
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FIG. 3: The map of the dynamics of distribution of a part 
of Gaussian wavepacket with a period of 130 fs and a 800 
nm wavelength in a photonic crystal oscillator with 100 pairs 
of layers and a 10% period variation. The spatial coordinate 
of system is taken over X axis, time is shown on Y axis and 
intensity of a signal is shown by color. 



Overall, this gives an equation 

M-l 

sinS M -i =0, £= Xi (9) 

The results of the calculation of the spectrum for ordi- 
nary one-dimensional photonic crystal and for one with a 
period variation are presented in Fig. 1,2. The localized 
levels between the upper bands of photonic bandgaps at 
the crystal edges and in the center of the structure are 
not exactly equidistant. This means that exact recon- 
struction of the form of the signal won't occur and the 
packet will irrevocably lose it's shape after some periods. 
However this fact is not crucial. Let us introduce the 



"non-equidistance parameter" rj = -^j? , where AK is an 
effective distance between levels and 5k is an average de- 
viation of the level interval from the AK. For an ideal 
oscillator -q = 0%, in case of ordinary photonic crystal 
rj = 10%, while in presence of adiabatic period variation 
rj = 2.5%. This means in a simple crystal the packet will 
lose it's shape approximately after 3 periods, while for 
the selected structure this will happen after 10 periods. 
This is enough for experimental purposes. 

One can consider a possible experimental realization 
concerning distribution of a wavepacket inside an ob- 
served structure. In common optical experiments, the 
sample is surrounded by air and placed on a substrate. So 
additional wide layers of air (n = 1) and silicon (n = 3.5) 
are added on the bounds of the system. We take a Gaus- 
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sian shaped wavepacket of a period of 130 fs and the 
wavelength 800 nm. These are common values for the 
state of the art Ti:Sap lasers. The optical lengths of 
the layers on the crystal edge are 4 = 200 nm. In this 
case the spectrum of the wavepacket lies inside the area 
of ~ 10 localized levels. The distribution of a signal is 
calculated in accordance of Eq.(p}. Fig. 3 shows the re- 
sults of calculation. After the initial moment, most of the 
pulse is reflected from the structure, because its carrier 
frequency corresponds to the photonic bandgap of crys- 
tal edge region. However, approximately 30% part of the 
signal in terms of intensity penetrates inside the oscillator 
due to the tunnelling and becomes localized. It starts a 
propagation with a periodical shape reconstruction. The 
localization due to the reflection of the wavepacket from 
the effective Bragg mirrors on the edges of the crystal 
leads to a periodical stops and changing of direction of 
the pulse propagation. The period of echo motion in thus 
created optical oscillator potential is 1 ps. 

Finally, let us discuss how the proposed optical echo 
effect can be observed in photonic crystal oscillator. One 



of the ways to detect the shape reconstruction is a time- 
resolved observation of the nonlinear response of the sys- 
tem. Indeed, near the stop-points of the oscillatoric mo- 
tion the pulse group velocity falls. Therefore, the pulse 
spatial length is smaller and its electromagnetic field is 
larger near these points, as can be observed from Fig. 3. 
Consequently, second or third optical harmonic response 
of the structure should look like a consequence of sharp 
peaks with a sub-picosecond period. This should be seen 
in the autocorrelation properties of the nonlinear signal. 
The third harmonic seems more perspective, as the effect 
in higher harmonics is more pronounced. On the other 
hand, there maybe an experimental problem concerning 
the absorbtion of the third-harmonic signal inside the Si- 
based structure at 800 nm pumpA^J^. This may require 
usage of the infrared pump laser (Ti-sapphire + para- 
metric oscillator system). 
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